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Abstract 

Pair production itself prevents the development of dyadospheres, hypothet- 
ical macroscopic regions where the electric field exceeds the critical Schwinger 
value. Pair production is a self-regulating process that would discharge a 
growing electric field, in the example of a hypothetical collapsing charged 
stellar core, before it reached 6% of the minimum dyadosphere value, keep- 
ing the pair production rate more than 26 orders of magnitude below the 
dyadosphere value. 
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1 Introduction 



Ruffini and his group [11121131111301110131101111111211131111113111 have proposed 
a model for gamma ray bursts that invokes a dyadosphere, a macroscopic region of 
spacetime with rapid Schwinger pair production [T2| , where the electric field exceeds 
the critical electric field value 

= — = ^ ^ 1.32 X 10^6 V/cm. (1) 
q nq 

(Here m and —q are the mass and charge of the electron, and I am using Planck 
units throughout.) The difficulty of producing these large electric fields is a problem 
with this model that has not been adequately addressed. Here I shall summarize 
calculations fH] showing that dyadospheres almost certainly don't develop astro- 
physically. 

The simplest reason for excluding dyadospheres is that if one had an astrophysical 
object of mass M, radius R > 2M, and excess positive charge Q in the form of 
protons of mass rUp and charge q at the surface, the electrostatic repulsion would 
overcome the gravitational attraction and eject the excess protons unless qQ < rripM 
or 

E qQ ^ TTlpM ^ Vflp ^ , o in-13 /^^O 



...p.,. ...p . 2 X in^^^ — (2) 
E, m2i?2 - ^2^2 ^ 4^2^ ^ • ^ \M J' ^ ^ 

where Mq is the solar mass. (If the excess charge were negative and in the form of 
electrons, the upper limit would be smaller by m/mp.) Then the pair production 
would be totally negligible. 

However, one might postulate the implausible scenario in which protons are 
bound to the object by nuclear forces, which in principle are strong enough to 
balance the electrostatic repulsion even for dyadosphere electric fields. Therefore, 
for the sake of argument, I did a calculation [TH] of what would happen under the 
highly idealized scenario in which the surface of a positively charged stellar core 
with initial charge Qo = M (the maximum allowed before the electrostatic repulsion 
would exceed the gravitational attraction on the entire core, not just on the excess 
protons on its surface) freely fell from rest at radial infinity along radial geodesies 
in the external Schwarzschild metric of mass M. 

This idealization ignores the facts that a realistic charged surface would (a) 
not fall from infinity, (b) have one component of outward acceleration, relative to 
free fall, from the pressure gradient at the surface, (c) have another component of 
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outward acceleration from the electrostatic repulsion, and (d) fall in slower in the 
Reissner-Nordstrom geometry if the gravitational effects of the electric field with 
Q ~ M were included. Because of each one of these effects, the actual surface would 
fall in slower at each radius and hence have more time for greater discharge than in 
the idealized model. Hence the idealized model gives a conservative upper limit on 
the charge and electric field at each radius, even under the implausible assumption 
that the excess protons are somehow sufficiently strongly bound to the surface that 
they are not electrostatically ejected. 

As we shall see, even in this highly idealized model, the self-regulation of the 
pair production process itself will discharge any growing electric field well before 
it reaches dyadosphere values. This occurs mainly because astrophysical length 
scales are much greater than the electron Compton wavelength, which is the scale 
at which the pair production becomes significant at the critical electric field value 
for a dyadosphere. Therefore, the electric field will discharge astrophysically even 
when the pair production rate is much lower than dyadosphere values. 

These calculations lead to the conclusion that it is likely impossible astrophysi- 
cally to achieve, over a macroscopic region, electric field values greater than a few 
percent of the minimum value for a dyadosphere, if that. The Schwinger pair pro- 
duction itself would then never exceed IQ-^^ t imes the minimum dyadosphere value. 

2 Schwinger discharge of an electric field 

In this section we shall analyze the pair production and discharge of an electric field 
produced by the collapse of the idealized hypothetical charged sphere or stellar core 
of mass M and initial positive charge Qo? assuming that somehow the excess charge 
on its surface is not electrostatically ejected, and assuming that the surface falls in 
as rapidly as possible, which is free fall from rest at infinity in an assumed external 
Schwarzschild geometry. 

As the surface radius R collapses, the electric field E = Qjr outside (r > R) 
produces pairs, with the positrons escaping and the electrons propagating in to the 
surface to reduce its charge Q. The pair production rate M per 4- volume [T7j is 
given by 




(3) 
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where 



w 



qE E qQ 
Anm'^MQ ( M \ (Qy^ f r ^ 



q x^y^Q, 



Mr:, KM J \2m) 



UM\fQV^lr_^^^ (4) 



B \MqI \MJ \2M 



B ^ -tAtt- ^ 42475. (5) 



f^yV M'<^5^^^^1.33xiO^ (6) 
Mp, \MJ \2MJ - Am?Mp, ^' 



A dyadosphere has E>Ec = m'^/q^w<7r 

IqJXm) \2M) 4m2A.Q 

For astrophysical electric fields anywhere near dyadosphere values, the elec- 
trons and positrons produced will quickly be accelerated to very near the speed 
of hght, so one will effectively get a null number fiux 4- vector n+ of highly rela- 
tivistic positrons moving radially outward and another null number fiux 4-vector 
n_ of highly relativistic electrons moving radially inward, with total current density 
4-vector j = qn^ — gn_. 

It is most convenient to describe this current in terms of radial null coordinates, 
say U and V, so that the approximately Schwarzschild metric outside the collapsing 
core may be written as 

ds^ = -e^^dUdV + r'^iU, V){de^ + sin^ Odcj)^). (7) 

Then Maxwell's equations (Gauss's law) gives 

47rj = 47rq{nldv - nidu) = V • F = ^ {Q^^'dy - Q'^'du) . (8) 

The 4-divergence of each of the number fiux 4- vectors n+ and n_ is equal to the 
pair production rate J\f, which leads to the following relativistic partial differential 
equation for the pair production and discharge process: 

Q,„. = -2.,.V"^=-^exp(-— , (9) 



or 



87r?rW = ^^^^ exp f-^) = ^(r'^^) = rU (rE) - r^EU Q , (10) 
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where is the covariant Laplacian in the 2- dimensional metric "^ds^ = —e^^dUdV 
and where □ is the covariant Laplacian in the full 4-dimensional metric. 

In my much more detailed paper I have analyzed this partial differential 
equation for the charge distribution over the entire spacetime region exterior to the 
charged surface and have found an approximation that reduces it to a relativistic 
ordinary differential equation for the evolution of the charge at the surface itself. 
However, here I shall confine myself to a Newtonian approximation to this evolution 
equation, which turns out to give results very close to the relativistic approximation. 

A convenient time parameter for describing the collapse of the surface is the 
velocity v the surface has in the frame of a static observer at fixed r when the 
surface radius R crosses that value of r. For free fall from rest at infinity in the 
external Schwarzschild metric of mass M, with r being the proper time along the 
surface worldlines, one gets 



dR 2M 

" = -d7 = V^' ^''^ 

so that the proper time remaining until the proper time at which the surface 
reaches the curvature singularity at i? = is — r = (4/3)M/f^ and R = 2M/v^ = 
(4.5M)^/^(rc — r)^/^. The Newtonian limit of this is when R ^ 2M, which gives 
f <^ 1 and T ^ t, the Schwarzschild coordinate time. In this limit, the surface moves 
negligibly during the time it takes for electrons to move inward from where they are 
created to the surface to reduce Q(t). 

li w = TiEc/ E is defined at each point outside the surface, let 

z(t) = wit, Rit)) = ' , = = — — — 12 

^ ' V , V E{R{t)) qQ{t) qQv^ ^ ' 

be the value of w at the surface itself. Since we shall find that the electric field E 
always stays far below the critical dyadosphere value E'c, we have z ^ 1, which will 
be used for various approximations below. 

Now the pair production rate (and assumed effectively instantaneous propagation 
of the electrons produced to the surface of the collapsing stellar core) gives 

-^^ -q U{r)ATxr^dr ^ -q — — / — e"^^ ^ --^-^^ — . (13) 
dt JR JR 2z-^e^ 

Then using R = 2M/v'^ and dv/dt ~ dr/dr = f^/(4M) leads to the following 

ordinary nonlinear first-order differential equation for z{v) in the Newtonian limit: 



V dz ^ 4 
z dv 
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-4 



1 - e 



-u 



(14) 
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^ -4:{z + 2) (l - e-^) + 5. {U 



where 

A = iM^I^a^ ~ 3.39643251 x lO^^, (15) 

71 

A-^. (16) 
U = z + \nz-\nA-2\nfi + 5\nv. (17) 

From Eq. (fT^ . one can see that the boundary conditions for Eq. fll4|) are that 
initially (r ^ t = — oo ■x^R = oo'^v = 0) v^z = 7rm^M^/(g(5) with the surface 
charge Q having its asymptotically constant initial value Qq, which will be taken to 
be its maximum allowed value, M, unless otherwise specified. One can see from this 
that both z and U start off initially at infinite values. The final value will be when 
the surface enters the event horizon of the black hole at R = 2M or f = 1. This 
is beyond the applicability of the Newtonian approximation being used here, but it 
turned out that the relativistic analysis |TH] gave very nearly the same answers. 

One can now differentiate Eq. (fTTj) for U{v), using Eq. (fT^. to obtain 

dU 
dv 

Since z ^ 1, this equation implies that U decreases to near zero (though it cannot 
reach zero, for if it could, the right hand side would be positive, contradicting the 
assumption that it dropped to zero and hence had a negative derivative on the left 
hand side). Then when [/ ^ 1, Eq. ()17p may be solved approximately for z{v) to 
give 

z{v) ^ InA + 21n/i — 51nf — 21n (In A + 21n/i — 51nw) 
^ ln/l-21nlnA+ (^1- (21n/i-51nt;) 

^ 57.33 + 1. 94 In /i + 4.85 In-. (19) 

V 

This then gives the ratio of the electric field at the surface to the critical electric 
field of a dyadosphere as being 

— = 1^ 0.0548 - 0.00185 In/x - 0.00463 In - < 1. (20) 

Ec z V 

One can improve this result by using a slightly improved formula (giving a 
roughly 2.5% correction for z ~ 57) for the radial integral in Eq. ()13|) for the 
pair production rate, by using a better explicit approximation for what U should 
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approach ai v = 1, and by using a numerical solution of the thus-corrected form of 
Eq. (fTTj) for z with this expression for U, to estimate that when /i = 1 and v = 1, 
z ^ 57.5843. Including the improved formula for the right hand side of Eq. (fT^ 
into the Newtonian differential equation (fT^ and then integrating it numerically 
from f = to = 1 for yU = 1 gave the result at the horizon of the solar black hole 
of 2; ~ 57.5845, so only the 6th digit changed from the algebraic estimate obtained 
without numerically solving the Newtonian differential equation. 

In ^H] I used a relativistic ordinary differential equation approximation to the 
partial differential equation and was able to deduce an explicit approximate 
relativistic result for fi = 1 and f = 1 of 2; ~ 57.60483, whereas the numerical 
solution of the relativistic ordinary differential equation gave z ~ 57.60480, differing 
from the explicit formula (using as input the values of m, q, and Mq) by only about 
one part in two million. However, I would estimate that the relativistic ordinary 
differential equation approximation to the partial differential equation would itself 
introduce absolute errors of the order of 10~^-10~^ in the value of z, so the numerical 
solution of that ordinary differential equation is not necessarily any better than the 
completely explicit approximate solution I also obtained. 

The difference between the Newtonian and the relativistic approximations for z 
on the horizon (f = 1) of a solar mass collapsing core (/i = M/Mq = 1) is about 
0.02, which I would guess is considerably larger than the error of my relativistic 
approximation (not given here, but in [18j), but it is still a relative difference of only 
about one part in three thousand for the idealized upper limit of the value of the 
electric field of a hypothetical charged stellar core collapsing into a black hole after 
falling in freely (no nongravitational forces on the surface) from starting at rest at 
radial infinity with the external metric being Schwarzschild. 

We can also give a heuristic derivation of the Newtonian result in the following 
way: We expect the self-regulation of the electric field to make z = nEc/E = 
•Kw^R^ / (qQ) change slowly, so 

1 dO 2 dR 1 /2M\ t , , 

(21) 



Q dt R dt M \ R J M' 

Then the pair production rate per 4- volume. A/" = (m'^/47r)e~'^/w^ ~ (m'^/47rz^)e~^'~^/^^ 
for z ^ 1 decreases roughly exponentially with r with e-folding length Ar ~ R/ {2z). 
Thus 

IdQ _ q , ^ 2j q 



r Af(r)A7Tr^dr ^ -^Af(r = RUnR^Ar 
JR Q 



Q dt Q Jr ^ ' Q 
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Equating this to -{1/ M){2M / Rfl"^ gives 

^ n ~ \2m) 



^5 ' 



(23) 



which imphes z + 2lnz ~ lny4 + 21nyU — 51nf, just as we got from the approximated 
solution of the ordinary differential equation in the Newtonian approximation. 

A dyadosphere would have E > E^, which implies z = tcEc/E < Zc = tt and 
^f = {m'^/4:7r)e~'/z^ >K = m^e'^ / {Atx^). But 
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,5 



7]-3g7r^5 



■Mr. 



0.672 X 10 



-26 



/i2 



(24) 



so the heuristic estimate gives the maximum pair production rate more than 26 
orders of magnitude below that of a dyadosphere. By comparison, the numerical 
solution of the approximate relativistic ordinary differential equation JHI gave, at 
v = h 



M 0.661168 X 10^26 

i 



0.0005545 In/i - 0.00001759 In^ /i). 



(25) 



about 1.6% less than the heuristic estimate gives. 

One can also calculate the maximum efficiency for converting the collapsing 
stellar core mass M into outgoing positron energy, 

QdQ 2/i2 ri 



MR 52 



' dv 
z^ 
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m 



lnA+-lnB - -lnln(AB^) + -Inu 
4 4 ^ ^ 4 

13^j3 



-1/2 



471QI 



3M\ q 



1.9 X 10" 



2 X 10" 



In 



2107^9^3 



In In 

4 



1 + 

m) 



In [M/Mq] 



3/2 



300 



-1/2 



'AT 



1/2 



m) 



3/2 



(26) 
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By comparison, the numerical solution of the approximate relativistic ordinary 
differential equation ^H] gave the more precise result 



The dominant factor in an estimate of the coefficient 0.0001855 (the upper limit on 
the efficiency if Qo = = ^0) is one-third the ratio of the electron mass to the 
proton mass, which is 0.0001815. This efficiency is too low for the pair production 
from these idealized collapsing charged cores to explain gamma ray bursts, even 
if it is admitted that this very implausible scenario (of the excess charge Q ~ M 
not getting expelled from the collapsing core by the huge electrostatic forces on it) 
comes nowhere near being able to form a dyadosphere. 

One can also calculate JHI that the probability of one of the particles annihilat- 
ing with an antiparticle is less than 10"^^, so the direct interactions of individual 
particles is negligible, consistent with what was assumed above. 

3 Conclusions 

If protons are bound to a collapsing stellar core purely gravitationally, the maximum 
electric field is more than 13 orders of magnitude below dyadosphere values: -Emax < 
1.2 X 10-^^{Mq/M)Ec. 

If protons are much more strongly bound, -Emax < 0.055£'c and A/'max < lO^^^A/'c, 
where Afc is the minimal dyadosphere production rate. 

The energy efficiency of the process for M ~ M© is very low, e ~ 1.86 x 



If one relaxed the assumptions of this model, such as the spherical symmetry, 
one would expect to get similar results, perhaps changing the pair production rates 
by factors of order unity that depend upon the precise geometry. However, it seems 
very unlikely that any modification could increase the maximum possible pair pro- 
duction rate by any significant fraction of the 26 orders of magnitude that the model 
above fails to achieve dyadosphere values. Therefore, the example analyzed strongly 
suggests that dyadospheres do not form astrophysically. 




(27) 



lO-^(M/M0)i/2(Qo/M)3/2 ^ (Qo/M)3/yM/(2.9 x lO^M©). 
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